Weak Ferromagnetic Exchange and Anomalous Specific Heat in ZnCu3(OH)6Cl2 
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Experimental evidence for a plethora of low energy spin excitations in the spin- 1/2 kagome antifer- 
romagnet ZnCu3(OH)6Cl2 may be understandable in terms of an extended Fermi surface of spinous 
coupled to a U(l) gauge field. We carry out variational calculations to examine the possibility 
that such a state may be energetically viable. A Gutzwiller-projected wavefunction reproduces the 
dimerization of a kagome strip found previously by DMRG. Application to the full kagome lattice 
shows that the inclusion of a small ferromagnetic next-nearest-neighbor interaction favors a ground 
state with a spinon Fermi surface. 
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The spin-1/2 Heisenberg antiferromagnet on a kagome 
lattice is a proving ground for the existence of a two- 
dimensional spin liquid. Can the combination of low 
spin, low coordination number, and geometric frustra- 
tion lead to spin disordered ground states and excita- 
tions with fractional statistics? If so, then it remains 
to establish whether spin-spin correlations are short- 
ranged with a gap to triplet excitations or whether the 
correlations are quasi-long-ranged with gapless excita- 
tions. Experiment is the ultimate arbiter of these fun- 
damental questions. Thus, the recent experimental re- 
alization of a spin-1/2 kagome antiferromagnet (KAF), 
ZnCu3(OH)6Cl2 also known as Herbertsmithite[T], is of 
great interest. Magnetic susceptibility, muon spin rota- 
tion, and NMR measurements show no evidence of mag- 
netic order [21 |3]. The bulk spin susceptibility exhibits 
Curie behavior down to temperature T — 0.5J, where 
J w 170 K, but then increases sharply with decreasing 
temperature, possibly saturating at T = 0. Likewise, the 
specific heat shows significant enhancement at low tem- 
peratures. At very low temperatures, the temperature 
dependence can be roughly fit to a power law with an 
exponent as small as a = 1/2. Over the range T = 5 
to 30 X 10""* J, a best fitg] to C„ cx T" yields a = 2/3. 
A small magnetic field suppresses the enhancement, sug- 
gesting that the contributing excitations are magnetic in 
origin. 

Various explanations for the experiments have been 
proposed. For instance, it has been observed that sam- 
ples of Herbertsmithite contain a significant number 
of impurties. Saturation of the susceptibility may be 
a Kondo effect due to impurity substitution of non- 
magnetic zinc ions by copper atoms [5] . A separate con- 
tribution from isolated spinful inpurities has been argued 
to explain the measured bulk susceptibilityfB]. Further- 
more, the temperature variation of the specific heat in an 
applied magnetic field resembles a Schotty anomaly and 
is similar to that attributed to defects in the spin-1 chains 
YzBaNiOs and NENP[I,i,. Powder NMR measurements 
of the local spin susceptibility suggest a second possibil- 
ity. Although the averaged local susceptibility tracks the 



bulk, regions of the NMR spectrum actually show a de- 
crease in the susceptibility with decreasing temperature. 
The discrepancy between local and bulk measurements 
may indicate the presence of impurities or may instead be 
attributed to the sampling of the susceptibility along dif- 
ferent and distinct crystallographic directions |3] [3] . If the 
latter case, the low temperature ferromagnetism may be 
due to a Dyzaloshinskii-Moriya (DM) anisotropv 1 Ol ITT] . 

To make a clear prediction for the specific heat at 
the lowest experimental temperatures, we use Gutzwiller- 
projected mean-field theory. Our approach is motivated 
in part by the appearance of a specific heat exponent 
a < 1 found for an extended Fermi surface of spinous 
coupled to a U(l) gauge field [T2]. A related composite- 
fermion theory has had some success in explaining the 
behavior of the half-filled Landau level 13]. Once the 
fermions are integrated out, the contribution of the gauge 
field to the specific heat can be obtained from the free 
energy [121 El!- In the absence of long-ranged density- 
density interactions between the fermions, it takes the 
form 

C(T) cx T^/^ (1) 

and thus, dominates conventional contributions such as 
that due to phonons. This signature provides motivation 
to investigate candidate states that possess large spinon 
Fermi surfaces. 

In the case of the nearest-neighbor (NN) Heisenberg 
model, our calculations reproduce the results of Ref. EJ 
namely, we find that the lowest energy state has no bro- 
ken symmetries and the corresponding mean-field state 
exhibits Dirac fermions at nodal points. Under the hy- 
pothesis that the experimentally observed weak ferro- 
magnetism is intrinsic to the pure lattice of Cu^"*" spins, 
we investigate the effect of adding a small next-nearest- 
neighbor (NNN) interaction with negative (ferromag- 
netic) coupling 

= Ji ^ S, • Sj - J2 ^ S, • S, . (2) 

For convenience, we set Ji = 1 in the following. We 



note that establishing the sign of the NNN spin-spin in- 
teraction from first principles requires a detailed under- 
standing of the quantum chemistry of ZnCu3(OH)6Cl2. 
An organic spin- 1/2 KAF has been argued to have mixed 
ferromagnetic NN and antiferromagnetic NNN bonds jTS] , 
but we leave open the possibility that other sub-leading 
terms, like the DM anisotropy, may also be important in 
a realistic model of Herbertsmithite. 

A rational approach to the systematic construction 
of variational spin liquid states can be realized by gen- 
eralizing the ordinary SU(2) Heisenberg antiferromag- 
net to SU(N) spins living in the self-conjugate repre- 
sentation, meaning that each site has N/2 fermions|16j. 
The resulting model may then be solved exactly in the 
large-N limit by decoupling the spin-spin interaction 
with the Hubbard-Stratonovich field Xij = {fl'^fja)- 
Appropriate variational states are then obtained by 
the usual Gutzwiller-projection procedure [T7]- To com- 
pare candidate variational wavefunctions in the physi- 
cal SU(2) limit, NN and NNN correlations are sampled 
by a standard determinantal wavefunction Monte Carlo 
algorithm [T51 [TH] that exactly enforces the single occu- 
pancy constraint //ViT + //Va = 1- The Xij field, 
interpreted as a bond amplitude order parameter, may 
then be adjusted to optimize the ground state energy to 
test for various dimerization instabilities. This approach 
has been shown to reproduce qualitatively the phase di- 
agrams of certain spin chains that are known exactly by 
other means [TT]. 

To test the reliability of this approach, we first study 
the nearest-neighbor Heisenberg model on a 2-leg kagome 
stip. Because of its one-dimensional character, the 
ground state has been accurately described by the es- 
sentially exact numerical density matrix renormalization 
group (DMRG) method where it was shown (20j to be 
strongly dimerized with five different values for the NN 
spin-spin correlations (S^ • Sj). The strip employed in 
the DMRG calculation has open boundaries that are 
terminated in such a way to remove triplet excitations 
from the low energy spectrum; nevertheless the nearest- 
neighbor spin-spin correlations are well converged in the 
interior. The NN spin-spin correlations (a) through (e), 
as labelled in Fig. [T] are -0.16429, -0.63506, -0.07087, 
—0.52862 and —0.33122 respectively HI], and the ground 
state energy per site based upon these interior values is 
Eq = -0.45929. 

For the variational calculation we adopt periodic 
boundary conditions. There are three distinct bonds 
per unit cell. The possibility of spontaneous symme- 
try breaking is then tested[17j by examining nontriv- 
ial patterns of gauge flux threading through triangles 
and hexagons. States that break time-reversal symme- 
try always appear to have high energy after projection; 
hence we follow Ref. and restrict attention only to 
fluxes and tt. NN correlations are sampled on a strip 
with length 100 sites and 250 total spins. We find that 
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FIG. 1: The kagome strip. Top: The optimal variational 
state. Numerical values indicate the relative hopping magni- 
tudes along the links and the value of the flux through each 
hexagon. There is no flux through the triangles. Bottom: The 
magnitude of the five symmetry-unrelated nearest-neighbor 
spin-spin correlation functions (S^ • Sj), labeled (a) through 
(e), is indicated by the thickness of the lines. See text for the 
actual values. 

among states that at most double the size of the unit 
cell, the lowest energy state has a flux of tt through ev- 
ery other hexagon and zero flux through the remaining 
hexagons and all triangles; thus, the variational state ex- 
hibits dimer order. The optimal arrangement of fluxes 
and bond magnitudes is depicted in the top panel of Fig. 
[T] Along the two legs, alternating bond magnitudes of 
approximately 1.1 and 0.9 yield the lowest energy, with 
the remaining third bond fixed at unit hopping ampli- 
tude. The doubled unit cell supports five symmetry- 
incquivalcnt values of the NN spin-spin correlation func- 
tion (Si • Sj) in agreement with the DMRG results. 
Along bonds (a) through (e) the spin-spin correlation 
functions are -0.08649(48), -0.67460(40), -0.06213(42), 
-0.46876(30), and -0.41077(36). The pattern shows 
good qualitative agreement with the DMRG calculation: 
the five distinct spin-spin correlation functions appear 
in the same order by size, and the variational energy 
per site of So = -0.45239(3) is only 1.5% higher than 
the DMRG value. Further relaxation of the energy to- 
wards the DMRG value may be expected if all five dis- 
tinct bonds were to be adjusted in magnitude. 

As the variational approach appears to give a good de- 
scription of the kagome strip, we now turn to the full 
2D kagome lattice. Several candidate spin liquid states 
have been examined [S . One competitive state has flux 
TT through the hexagons and no flux through the trian- 
gles. This TT-hexagon state has Dirac nodes at the Fermi 
level[5 . Another state of interest has no flux penetrat- 
ing through the lattice anywhere. As shown in Fig. [2] 
this uniform state possesses an extended Fermi surface 
at half-filling. 

The energies of the candidate states are calculated on 
an oblique lattice with 12 x 12 3-site unit cells (432 total 
sites). On finite lattices, the choice of boundary con- 
ditions is important. In the case of periodic boundaries. 
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FIG. 2: (Color online) Mean-field spinon dispersion of the 
uniform phase of zero flux, (a) Complete dispersion showing 
that one of the three bands is completely flat, (b) Middle 
band cut by a plane at the Fermi energy corresponding to 
half-filling. Note that the Dirac nodes are not located at 
the Fermi energy, (c) Rotation of (b) to highlight the nearly 
circular Fermi surface. 



FIG. 3: Two possible dimer orderings with unit cells outlined 
by dashed parallelograms. Thick lines indicate larger hopping 
magnitudes and hence, larger NN spin-spin correlations. Top 
panel: Dimer order that maximizes the number of hexagons 
with exactly three dimers [24]. Bottom panel: Plaquette order 
that outlines the shape of a star|22j. 



there is an ambiguity in the fiUing of independent-particle 
states at the Fermi level, as only some of the degener- 
ate independent-particle states are filled. To resolve the 
ambiguity, the degeneracy must be lifted. Mixed peri- 
odic and antiperiodic boundaries accomplishes this [53] 
at the cost of breaking lattice rotational symmetry. Con- 
sequently, the mixed boundary conditions induce modu- 
lation in observables such as (Sj • S^). The modulation, 
however, decreases with increasing system size, vanish- 
ing in the thermodynamic limit. For the relatively large 
lattice studied here, the modulation is smaller than the 
uncertainty in the calculated spin-spin correlation func- 
tion due to the statistics of finite samples. 

We check for instabilities of the variational states to- 
ward two different types of dimer patterns. The first 
pattern is a dimer order of the type first proposed in 
Ref. l24l based on 1/7V corrections to the large-N solution. 
Dimer coverings that maximize the number of hexagons 
with exactly three dimers ( "perfect hexagons" ) are ener- 
getically preferred because the perfect hexagons resonate 
like the alternating single and double bonds of a ben- 
zene ring. As the resonance is a second-order process, 
it lowers the energy. There can be at most one perfect 
hexagon for every 18 sites[21]. The top panel in Fig. 
[3] shows an 18-site unit cell with one perfect hexagon. 
Larger unit cells support other possible dimer coverings 
that maximize the density of perfect hexagons including 
the honeycomb pattern shown in Fig. 1(b) of Ref. |24] 
and the striped pattern shown in Fig. 5(b) of Ref. [55] 
For the NN Heisenberg model, a systematic dimer ex- 
pansion determines that the honeycomb and stripe pat- 



TABLE I: The lattice-averaged spin-spin correlation function 
{Si ■ Sj) for an oblique lattice of 12 x 12 3-site unit cells (432 
sites total) with mixed periodic and antiperiodic boundary 
conditions. The magnitudes of the bond amplitudes are mod- 
ulated by 5% for the ordered states. Numbers in parentheses 
are the statistical errors due to the finite number of samples. 



TT-hexagon state 
dimer order 
plaquette order 

uniform state 
dimer order 
plaquette order 



e NN 

-0.42866(2) 
-0.42856(3) 
-0.42841(4) 
-0.41215(2) 
-0.41205(2) 
-0.41203(2) 



{i,j) e NNN 

-0.02297(21) 

-0.02277(22) 

-0.02271(17) 

0.07724(13) 

0.07873(10) 

0.07701(13) 



terns have ground state energies per site of -0.432088216 
and -0.4315321, respectively [55]. These energies compare 
quite well to the ground state energy of a 36 site clus- 
ter, -0.438377, found by exact diagonalization j27j . The 
second dimer configuration has plaquette order of a type 
that has been proposed as a potential instability of the 
TT-hexagon state. As shown in the bottom panel of Fig. 
[3] the plaquette takes the shape of a star formed by the 
outer edges of six triangles that share a common hexagon. 
This pattern of bonds would be induced by the dynam- 
ical formation of a non-chiral mass term that gaps out 
the nodal fermions |22). 

Table [T] summarizes the variational results. Instabili- 
ties are tested by imposition of a small (5%) modulation 
in the magnitude of the Xij bond amplitudes. In addi- 
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tion to the tabulated results, we also find that imposi- 
tion of a 5% modulation of the honeycomb type upon 
the TT-hexagon state raises the energy only slightly, to - 
0.42860(3), while at the same time gapping out the Dirac 
nodes. This result lends support to the idea[Ml 
that the nearest-neighbor antiferromagnet is actually 
dimer ordered. Furthermore, all of the dimer cover- 
ings we tested that maximize the number of perfect 
hexagons are energetically preferable to the star plaque- 
tte pattern [5J 122]. We speculate that further variational 
tuning will yield an ordered state of perfect hexagons 
that is energetically competitive with the best estimates 
for the ground state energy of the NN spin-1/2 KAF. 

Because the NNN spin-spin correlations are negative in 
the TT-hexagon phase but positive in the uniform phase, 
the uniform state is energetically favored over the tt- 
hexagon state for sufficiently large ferromagnetic NNN 
coupling, J2 > 0.16 (the ground state remains antiferro- 
magnetic up to much larger values of J2). At J2 = 0.16, 
dimer modulation of the uniform state further lowers the 
energy, with a minimum in the energy occurring at ap- 
proximately 4% bond modulation. This 4% bond modu- 
lation induces a 18% modulation in the value of (Si • Sj) 
on symmetry-distinct bonds. Somewhat surprisingly, the 
mean-field theory continues to support gapless excita- 
tions despite the broken translational and rotational sym- 
metries, with little change in the density of states at the 
Fermi energy. This is because ordering wavevectors cor- 
responding to large (multiple of 18 sites) unit cells do not 
span the spin on Fermi surface. Thus, a small ferromag- 
netic NNN interaction favors a spin-liquid phase with a 
large spinon Fermi surface, suggesting that spinous can 
remain deconfined even in the presence of dimer order. 
We are not aware of any rigorous arguments (eg. Ref. 
that preclude this possibility. 

The large spinon Fermi surface exhibits at low temper- 
ature both a relatively large (but finite) magnetic suscep- 
tibility and an enhanced specific heat, Eq. [l} due to inter- 
actions between spinous and the U(l) gauge field. These 
low energy properties of the uniform state agree qualita- 
tively with those from available Herbertsmithite exper- 
iments, though to reproduce quantitatively the anoma- 
lously large susceptibility, it may be necessary to invoke 
additional contributions from impurities [6 and/or the 
DM interaction [inj E] . 
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